Observations of two-fold shell filling and Kondo effect in a graphene nano-ribbon 

quantum dot device 
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A graphene nanoribbon (GNR) with orientation along its principle axis was obtained through 
a mechanical tearing process, and a quantum dot device was fabricated from the GNR. We have 
studied the transport property of the GNR quantum dot device down to dilution refrigerator tem- 
peratures. Two-fold charging periodicity was observed in the Coulomb-blockade measurement, sig- 
naling a shell-filling process with broken valley degeneracy. In one of the smaller Coulomb diamonds, 
Kondo-Iike resonance were observed, with two conductance peaks displaced symmetrically from the 
zero bias voltage. The splitting of Kondo resonance at zero magnetic field suggests spin-polarization 
of the quantum dot, possibly due to the edge states of a zigzag GNR. 

PACS numbers: 72.80.Vp, 73.23.Hk, 73.22.Pr, 73.63.Kv 



The discovery of quantum Hall effect in graphene [l|, Q 
has triggered a lot of studies on this material. Atten- 
tion has also been paid to the nanoribbon forms of the 
material, due to their intriguing physical properties for 
basic research and potential applications in future nano- 

electronic circuits S li I^IM III, Isl, S, [id, llll llJ, E, [3, 

As an unrolled version of single- walled carbon nanotubes, 
GNRs could be flexibly tailored to various shapes with 
controlled width, length, orientation and even intra- 
molecular junctions that are suitable as the building 
blocks for electronic circuits. With the unrolled struc- 
tures, however, the edges of GNRs are expected to play 
an important role, introducing new physics different from 
carbon nanotubes. Theoretical studies 0, 0, H, H, 0, Q 
show that GNRs with zigzag edges have spin polarized 
edge states, and are half-metals under a large trans- 
verse electric field. Similar prediction has also been 
made for bilayer zigzag GNRs [ll[i3l- Experimentally, 
GNRs have been fabricated through a reactive ion etch- 
ing method with the help of electron beam lithography 
[121, or via a chemical method 



14l |. Electron trans- 



port measurements reveal that these GNRs are all semi- 
conducting, with a gap inversely proportional to their 



widthjl2l. |14|. In the experiments of Ref. [lj|, sub-10- 



nni wide GNRs along principle axes were also reported. 
However, there is no systematic experimental study on 
GNRs with identified edge orientations. 

In this paper, we report our electron transport mea- 
surement on a quantum dot (QD) device made of a very- 
long GNR along one of the principle axes of a graphene 
sheet. Coulomb blockade diamonds with two- fold charg- 
ing periodicity were observed. Moreover, inside one of the 
small diamonds a fine peak-dip-peak conductance struc- 
ture was discerned around zero bias voltage. We believe 
that these results indicate the breakdown of both the 
valley degeneracy and the spin degeneracy at different 
energy scales in the GNR. 



We use the "Scotch Tape" technique [18| to exfo- 



liate highly ordered pyrolytic graphite (HOPG) onto 
degenerate-doped p-type Si substrates with a 100 nm 
thick Si02 layer. Besides obtaining large pieces (several 
/im^ and up) of graphene occasionally we can get 
very long GNRs. Figure 1(a) is a scanning electron mi- 
croscope image of the GNR used in this experiment. It is 
about 12 microns long and seems to have a uniform width 
over its entire length. Atomic force microscope (AFM) 
shows that the GNR has a width of w = 60 nm and a 
height of ~ 1.3 nm. The height of the GNR indicates 
that it has at most two graphite layers [l^ . 

The fact that some of the mechanically exfoliated 
GNRs have a constant width over a very long distance 
would suggest that they have two extremely parallel 
and smooth edges. This is most possible when both 
edges are along one of the principle axes of the graphene 
sheet. The existence of favorable tearing directions was 
reported by Geim's group on large pieces of graphene [20| . 
and by Dai's group on some of their chemically derived 
GNRs[14]. As shown in Fig. 1(a), besides the very long 
GNR, there is another GNR whose direction changes ex- 
actly by 30 degrees across a "kink" structure, presumably 
from zigzag to armchair or vice versa. The "kink" pro- 
vides further evidence that mechanically exfoliated GNRs 
are likely along the principle axes. 

Pd electrodes of 50-nm thick, 500-nm wide and 500- 
nm apart were patterned onto the GNR using standard 
electron beam lithography techniques (Fig. 1(c)). Two- 
probe differential conductance between neighboring elec- 
trodes of the device was measured using low-frequency 
lock-in techniques with an ac modulation voltage of 10 
/iV. The measurements were performed in a Quantum 
Design PPMS system with temperatures T down to 2 K, 
and in a dilution refrigerator with T down to 20 mK. 

Figure 2(a) shows the differential conductance G = 
d//(i Vbias in a color scale as a function of bias voltage 
Vbias and back-gate voltage V^, measured from one of the 
segments of the GNR at T = 2 K in the PPMS. The 
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FIG. 1: (color online) (a) Scanning electron microscope im- 
age of mechanically-torn graphene nanoribbons (GNRs). The 
short GNR kinks exactly by 30 degrees, (b) Atomic force mi- 
croscope image of the long GNR near an electrode, indicating 
the GNR width of about 60 nm. (c) Optical image of the 
electrodes fabricated onto the long GNR (a line is drawn to 
represent the GNR). 



results measured from other segments are qualitatively 
the same, and are similar to those reported previously 

13, EH, 



22l | . The conductance is suppressed at small Vbias 
due to the Coulomb blockade (CB) effect, and the sup- 
pression is enhanced in the vicinity of the Dirac point 
centered at Vg ~3 V, where the carrier concentration is 
the lowest. In Fig. 2(b) we show the zero bias CB oscil- 
lations near the Dirac point in the range Vg = 3.24 — 3.46 
V, measured in a dilution refrigerator at T = 600 mK 
and 20 mK, respectively. The conductance peaks become 
much sharper at 20 mK. Figure 2(c) shows the stability 
diagram at this temperature, in the same Vg window as 
in Fig. 2(b). The most salient feature is the two- fold 
charging periodicity with varying Vg . The bigger CB di- 
amonds are about 23 meV wide, and the smaller ones 
are about 15 meV wide. The height of the smaller di- 
amonds is around 0.6 meV. Excited states at an energy 
scale of ~ 0.3 meV can be recognized at the edges of the 
diamonds, as marked by the arrows in Fig. 2(c). 

The appearance of CB diamonds with excited states in- 
dicates that the GNR forms a QD at low temperatures. 
The effective size of the QD can be estimated from the 
width of the smaller CB diamonds, AVg_s — e/Cg ^ 15 
meV, where e is the electron charge and Cg the capac- 
itance between the QD and the back gate. Assuming 
a simple parallel plate capacitance, the length of the 
GNR QD estimated is about 510 nm that happens to be 
very close to the distance between the contact electrodes 
L — 500 nm. However, a more accurate model taking into 
account the finite size effect of the GNR plate ^] gives an 
effective length L* « 175 nm. The difference between L 
and L* may reflect the existence of long depletion regions 
that serve as tunnel barriers near the electrodes 2J] . We 
rule out the possibility that the device contains multiple 
QDs, because no joint CB diamonds are observed. 

Apparently, the two-fold periodicity observed in Fig. 
2(c) resembles the shell- filling process of a bar-like QD, 
where the orbital degeneracy is lifted by the spatial con- 
finement while the spin degeneracy is still held. Shell- 
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FIG. 2: (color online) (a) Color image of differential conduc- 
tance, G = d//dVbias, as a function of bias voltage Vbias and 
back gate voltage V^, measured at T =2 K. (b) Coulomb oscil- 
lation of G at T = 600 mK and 20 mK. (c) Stability diagram 
of the GNR device at 20 mK and in a Ig window from 3.26 to 
3.46 V. Coulomb diamonds with a two-fold periodicity along 
with the excited states outside the diamonds (indicated by 
the arrows) are observed. 



filling feature is widely observed in QD devices made 
of carbon nanotubes 25, 2^, 27 1 or semiconductors [1^. 
When the carrier occupancy of the QD, N, is odd, addi- 
tion of one carrier only needs to overcome the Coulomb 
energy, whereas it has to pay an extra energy to occupy 
the upper shell when N is even. As a result, the size 
of the CB diamonds shows two-fold periodicity, with the 
small diamonds correspond to odd N and the larger ones 
correspond to even N. As it is known, the electrons 
in 2D graphene exhibit a four-fold degeneracy near the 
Dirac point, i.e., two-fold spin degeneracy plus two-fold 
K and K' valley degeneracy. The observation of a two- 
fold charging periodicity possibly signals the lifting of 
the valleydegeneracy due to the existence of edges in the 
GNR 1,1, H 

Another interesting feature of the data is the appear- 
ance of a Kondo-like conductance enhancement in one 
of the Coulomb diamond centered at Vg =3.35 V (the 
sixth diamond from the right in Fig. 2(c)). Since this 
diamond is smaller than its neighbors, the electron occu- 
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FIG. 3: (color online) (a) Color-scale plots of the conductance 
G in one of the Coulomb diamonds with odd number of carrier 
occupation (the sixth diamond from the right in Fig. 2(c)). 
Conductance peak-dip-peak structure near zero bias voltage 
and its evolution with Vg can be resolved. T = 20 mK. (b) 
Peak position as a function of Ig. (c) G — Vbias curves mea- 
sured at different temperatures. The peak-dip-peak structure 
disappears above ~ 400 mK. (d) The amplitude of conduc- 
tance enhancement AG in (c) as a function of temperature. 
Insets of (d) depict two different ways of estimating AG. 



pancy N on the QD is presumably odd, leaving an un- 
compensated electron spin on the QD. It seems that this 
uncompensated electron spin mediates a resonant screen- 
ing process between the conduction electrons in the elec- 
trodes, thus provides a conduction mechanism in the CB 
regime 



M Im lU, HI , as shown in Fig. 3. 



Kondo resonance usually appears as a single conduc- 
tance peak at zero bias voltage (Fig. 4(a)). However, the 
conductance enhancement observed here exhibits a peak- 
dip-peak structure around zero bias. The two peaks are 

0.1 meV apart at the most left side of the diamond, 
getting closer with increasing Vg, and apparently merging 
together in the right half of the diamond [s^. 

Usually the height of a Kondo conductance peak fol- 
lows a logarithmic temperature dependence below the 
Kondo temperature, and saturates gradually at very low 
temperatures. The temperature dependence data of the 
peak-dip-peak structure are shown In Fig. 3 (c). As can 
be seen in Fig. 3(d), the amplitude of the overall struc- 
ture, AG, extracted using two different ways (depicted 
in the insets), does roughly follow a logarithmic temper- 
ature dependence below 400 mK. 

It is known that the Kondo resonance peak splits into 
two peaks that are symmetrically displaced from the zero 
bias if there is a preferred spin orientation (i.e., when the 
spin degeneracy is lifted, as illustrated in Fig. 4(b)) [3l|. 
Such splitting has been observed in both 2D semiconduc- 



tor QDs and carbon nanotubes, where the spin degener- 
acy is lifted by either an applied magnetic field [3J|or 
by the exchange coupling to ferromagnetic contacts [35| . 
The lifting of spin degeneracy in this experiment is sur- 
prising in that there is neither external magnetic field nor 
ferromagnetic contacts. 

Although the exact nature is obscure, the observed 
peak-dip-peak structure may be qualitatively understood 
by assuming some degree of spin polarization on the GNR 
QD itself. As aforementioned, preferred spin orientation 
is expected at the edges of a zigzag GNR 0,[i[il,[l3- In 
following paragraphs, by analyzing the CB data in Fig. 2 
we will show that the GNR measured in this experiment 
is quite likely to be a zigzag one. 

Let us first show that the QD is unlikely made of an 
armchair GNR. The band dispersion of an armchair GNR 
takes the same Fermi velocity {vp — 1.0 x 10^ ni/s) as 
in a bulk graphene, except a gap near the neutral point 
in the semiconducting cases. Taking L* « 175 nm, the 
level spacing due to longitudinal confinement for this bar- 
like GNR is eL — hvp/'^L* « 12 meV, and the energy 
spacing between subbands is e^, = hvp/2w ~ 34 meV. 
In order to explain the excited state energy of A « 0.3 
meV as observed in the CB measurement, one has to 
assume that there are many subbands across the Fermi 
level. Simple estimation [3^ shows that the Fermi energy 
Ei? ~ 430 meV, and the number of carriers in the GNR 
N ~ 710. This would lead to a Dirac point (the neutral 
point) that is located 20 V away from the Vg window of 
Fig. 2 (c). However, the experimental data indicate that 
the Dirac point is located at several volts of V^, or at 
worst within 10 V from the Vg window. In addition, if 
the GNR is armchair, in the regime near the Dirac point 
there should be at least a few big CB diamonds whose 
height is comparable to the confinement energy spacing 
of 12 meV and whose width is wider than 250 meV in Vg. 
However, as shown in Fig. 2 (a), there is no CB diamond 
having these large sizes. 

On the other hand, for a zigzag GNR, the energy 
bands near the Fermi level is nearly flat due to the ex- 
istence of edge states. Therefore, the confinement along 
its length direction yields a relatively small energy level 
spacing, so that the lack of large CB diamonds is nat- 
urally understood. To estimate the level spacing (the 
excited states energy) , we recall that the bands of edge 
states are slightly dispersive because of the next nearest- 
neighbor interaction , with a small group velocity Vp 
at the order v'p/vp ~ t'/t ~ 0.03, where t « 3.0 eV and 
t' 0.1 eV are the first and the second nearest-neighbor 
hopping energy [ill, respectively. The level spacing es- 
timated from the dispersion of edge states is roughly 
e'p « hv'p/2L* ~ 0.4 meV, consistent with the observed 
values. Therefore, we believe that our GNR is most likely 
to be a zigzag GNR rather than an armchair one. 

Based on the spin configuration of the edge states in a 
zigzag GNR 0, Q that the spins on the same edges are 
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FIG. 4: (color online) (a) Elastic Kondo resonance in a quan- 
tum dot device provides a conduction mechanism at zero bias 
voltage in the Coulomb blockade regime, (b) Inelastic Kondo 
resonance due to lifted spin degeneracy leads to a pair of 
conductance peaks at finite bias voltages (symmetrical about 
zero bias), (c) A toy model depicts the possible Kondo res- 
onance in a zigzag GNR. The spins are parallel on one edge 
but are antiparallel between opposite edges. Kondo resonance 
becomes inelastic when the electrostatic environment of the 
two edges is different. 



parallel but on opposite edges are antiparallel, we present 
a toy-model in Fig. 4 (c) to account for the lifting of spin 
degeneracy. In this model, electrons with different spin 
will fill onto different edges. When the local electrostatic 
environment of the two edges are different, filling to dif- 
ferent edge will cost different energy. This would lead to 
the observed spin splitting in odd N GNR QD. In this 
scenario, the gate- voltage dependence of the Kondo peak 
splitting can also be naturally explained through the de- 
pendence of local electrostatic environment on the gate 
voltage. In fact, it is almost unavoidable to have ran- 
domly trapped charges in Si02 surrounding the GNR. 
The distribution of these charges, thus their stray field, 
is tunable by applying a gate voltage. To further clar- 
ify the issue, more experiments such as magnetic field 
dependence of the peak splitting arc demanded. 

To summarize, we have observed two-fold charging 
periodicity in a QD device made of a GNR. We also 
observed Kondo-likc conductance peaks in one of the 
Coulomb diamonds. The observation of Kondo splitting 
at zero magnetic field suggests the lifting of spin degener- 
acy due to spin polarization on the QD itself, in favor of 
theoretical predictions on a zigzag GNR where the spin 
polarization is due to its highly degenerate edge states. 
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